Cosmological Tests for a Linear Coasting Cosmology 

Abha Dev* Margarita Safonovaj Deepak Jainf and Daksh Lohiya^ 
Department of Physics and Astrophysics, 
University of Delhi, Delhi-110007, India 

February 5, 2008 

(N : 

o ■ 
o ■ 

CN . Abstract 
C ' 

^ ' A strictly linear evolution of the scale factor is a characteristic feature in several 

■ classes of alternative gravity theories. In this article we investigate the overall viability 
! of an open linear coasting cosmological model. We report that this model is consistent 

with gravitational lensing statistics (within lex) and accomodates old high-redshift 
J> ■ galaxies. We finally conclude that such a linear coasting, a{t) = t, is not ruled out on 

O , basis of these observational tests. 

in . 
^ ■ 

O ; 1 Introduction 
(N . 

■ Standard cold dark matter (CDM) cosmology presents serious theoretical and observational 
P_|. problems as a model for an acceptable description of the Universe. This has motivated a 
Q ! search for alternative cosmological models [0. The first problem comes from the conflict 
^ I between age of the Universe and the age of the oldest stars in Galactic globular clusters. 
^ ■ The age constraints from old galaxies at high redshifts render "the age problem" even more 
^ ■ acute 0. 

"-j , Another major difficulty is the cosmological constant (A) problem. While it is difficult to 

', have a decent theoretical justification for A [Q, there are serious doubts on the compatibility 
I of the Hubble Deep Field (HDF) data with large A. Initial analysis of Maoz and Rix ^ 
placed a somewhat stringent upper limit of Qa < 0.7 which have been somewhat softened 
by later analysis of Cooray et al. p. 

These problems have generated a lot of interest in an open FRW model with a linear 
evolution of the scale factor, a{t) cx t. In such a cosmology the universe expands with 
constant velocity; hence the term coasting cosmology ||^. Notable among such models is a 
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recent idea of Allen in which such a scaling results in an SU{2) cosmological instanton 
dominated universe. The Weyl gravity theory of Manheim and Kazanas 0] makes space for 
yet another possibility. Here again the FRW scale factor approaches a linear evolution at 
late times. 

The need for investigating such a model comes from several considerations. Particle 
horizons occur in models with a{t) ^ for a < 1 and, therefore, linear coasting model 
does not suffer the horizon problem. Also, linear evolution of a scale factor is supported in 
alternative gravity theories (eg. non- minimally coupled scalar-tensor theories), where it turns 
out to be independent of the matter equation of state (see p and references therein). The 
scale factor in such theories does not constrain the matter density parameter and, therefore, 
does not present any flatness problem. Moreover, the age of the coasting universe is 1.5 times 
the age of standard CDM universe 0, which makes this model comfortably concordant 
with the ages of globular clusters. Finally, a linear coasting cosmology, independent of 
the equation of state of matter, is a generic feature in a class of models that attempt to 
dynamically solve the A problem ||^. Such models have a scalar field non- minimally coupled 
to the curvature of the universe. With the evolution of time, the non-minimal coupling 
diverges, the scale factor quickly approaches linearity and the non-minimally coupled field 
acquires a stress energy that cancels the vacuum energy in the theory. 

Interestingly it was noted by Perlmutter et al. that the curve for = = (for 
which the scale factor would have linear evolution) is "practically identical to the best fit 



plot for an unconstrained cosmology" . Recently, it was shown by Dev et al. [|TT| that open 
linear coasting cosmology presents a good fit to the SNe la data. It was also demonstrated 
that this model is consistent with the primordial nucleosynthesis |0 . 

In this paper we consider constraints on the index a of a power law cosmology, a{t) cx t", 
from two different tests: gravitational lensing statistics and age estimates of old high-redshift 
galaxies. The expected frequency of multiple imaging lensing events is a sensitive probe for 
the viability of a given cosmology. In view of the successful results of the above mentioned 
works ||1T], [12|, we used this test to constrain the power index a of the scale factor. Expected 
number of lens systems depends upon the index a through the angular diameter distances. 
By varying a, the number of lenses changes which on comparison with the observations gives 
us the constrain on a. Age measurements of old high-redshift galaxies give lower bound on 
the power index a. This is based on the fact that the age of the universe in a given redshift 
is greater than or at least equal to the age of its objects. Since the age of the universe is a 
function of a, we find the value of a which permits the existence of these old galaxies. 

In Section 2 we introduce the ansatz for the power law cosmology and derive the angular 
diameter distance formula. In Section 3 we describe the tests and the constraints they present 
on the power index for that cosmology. We summarize our results and present the result of 
a joint constraint from various observational tests in Section 4. 
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2 Linear Coasting Cosmology 



We consider a general power law cosmology with the scale factor given in terms of two 
arbitrary dimensionless parameters B and a 

t 



a{t) = B 



Ho \to 



for an open FRW metric 



dr"^ 



1 + r' 



1 



(2) 



Here t is cosmic proper time and r, 6, (p are comoving spherical coordinates. 

The expansion rate of the universe is described by a Hubble parameter, H{t) = a/ a = a/t. 
The present expansion rate of the universe is defined by a Hubble constant, equal in our model 
to Hq = a /to (here and subsequently the subscript on a parameter refers to its present 
value). The scale factor and the redshift are related to their present values by a/ao = {t/to)". 
As usual, the ratio of the scale factor at the emission and absorption of a null ray determines 
the cosmological redshift z by 

(7,(1 

l + z , 



a{z) 



and the age of the universe is 



t{z) 



a 



Hoil + zYl- ■ 



Using 



), we define the dimensionless Hubble parameter 

Hiz 
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The present 'radius' of the universe is defined as (see Eq. 

c 



ao = B 
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In terms of the parameters a and B, the angular diameter distance between two different 
redshifts is: 



DA(zi,Z2,a) 
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In a limiting case, a — > 1, we obtain the following expression 



BAzx-.zi) 
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[l + z^ 
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;i + zi)VB(i + z2)^ 



The look-back time, which is the difference between the age of the universe when a 
particular light ray was emitted and the age of the universe now, we find as 

edt c 

dzL Hoil + z-L)"^ ' 



(9) 
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3 Testing the model against observations 



3.1 Gravitational lensing statistics 



We consider a sample of 867 {z > 1) high luminosity optical quasars which include 5 lensed 
quasars (1208 + 1011, H 1413 + 117, LBQS 1009 + 0252, PG 1115 + 080, 0142 + 100). 
This sample is taken from optical lens surveys such as the HST Snapshot survey the 
Crampton survey , the Yee survey |]15| , Surdej survey |jl6[ , the NOT Survey [0 and the 
FKS survey ||T^. The lens surveys and quasar catalogs usually use V magnitudes, so we 
transform my to a B-band magnitude using an average B — V colour of 0.2 mag [|19|, |2^, ^ 
The differential probability dr of a beam having a lensing event in traversing dzi, is 



dr = noil + zlY a ^dzi^ . 

dz-L 



(10) 



where no is the present comoving number density of the lenses, a is the cross-section for 
lensing event and the quantity cdt/dzi^ is given by (|^). 

For simplicity we use the Singular Isothermal Sphere (SIS) model for the lens mass 
distribution. The cross-section for lensing events for the SIS model is given by |2^ 



a 



16 - 



OS 



(11) 



where v is the velocity dispersion of the dark halo of the lensing galaxy. We define -Dql, -Dqs 
and Dls as the angular diameter distances from the observer to the lens, to the source and 
between the lens and the source, respectively. 

Assuming no evolution of the galaxies, the comoving number density is modeled by 
Schechter function as 



$(L, z = 0)dL 



L\ dL 



(12) 



where 0,,, a and are the normalization factor, the index of faint-end slope, and the 
characteristic luminosity, respectively. We also assume that the velocity dispersion of dark 
matter halo v is related to the luminosity L by the Faber- Jackson relation for E/SO galaxies 
v = v, [LjUfl^ 

The differential probability is given byp3 



dr 



F*(1 + zl)'— . 

c V 



Ho f HoDolDi^\^ cdt_ 
cDqs ) dzi^ 
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where 
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(14) 



is the dimensionless quantity, which measures the effectiveness of matter in producing mul- 
tiple images. Table 1 lists Schechter and lens parameters for E/SO galaxies, as suggested by 
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(hereafter LP EM parameters). We neglect the contribution of spirals as lenses, as their 
velocity dispersion is small in comparison to E/SO galaxies. 

The constraints obtained on the cosmological paramters are highly dependent on the 
choice of lens and Schechter paramters. The lens and Schechter parameters should be deter- 
mined in a highly correlated manner from a galaxy survey. The use of parameters derived 
from various surveys might introduce error. We consider LPEM parameters as they form one 
such set of paramters and they also take into account the morphological distribution of the 



E/SO galaxies 1 25 1 . Recently several galaxy surveys have come up with a much larger sample 



of galaxies. This has improved our knowledge of the galaxy luminosity function. However, 



these surveys donot classify the galaxies by their morphological type ||2^, |2^, ^ . 

The differential optical depth of lensing in traversing dzi^ with angular separation between 
and (j) + d(f) reads as 



dzi^dcj) 
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where = A6'/87r(f^=/c)^ and the velocity dispersion corresponding to the characteristic 
luminosity in ([12|) . 

We make two corrections to the optical depth to get the lensing probability: magnification 
bias and selection function. Magnification bias, B(m, z), is to take into account the increase 
in the apparent brightness of a quasar due to lensing, which, in turn, increases the expected 
number of lenses in flux limited sample. 

The bias factor for a quasar at redshift z with apparent magnitude m is given by [^, ^ 
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B(m, z) = Ml B{m, z, Mq, M^) , (16) 
where 

fdNoY^ r^^dMdNn, . . . 

B(,n, z, M.) = 2 (^] (,„ + 2,5 log(M). .) . (17) 

In the above equation {dNQ{m, z)/dm) is the measure of number of quasars with magnitudes 
in the interval (m, m + dm) at redshift z. We can allow the upper magnification cutoff M2 to 
be infinite, though in practice we set it to be M2 = 10^. Mq is the minimum magnification 
of a multiply imaged source and for the SIS model Mq = 2. 

We use Kochanek's "best model" (K96) for the quasar luminosity function: 



dN. 



[m, z) oc (^io-"(") + io-^('"-™)) ^ , (1^ 



dm 

where the bright-end slope a and faint-end slope h are constants, and the break magnitude 
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m evolves with redshift: 



m 



rrio + (-2 — 1) for z < 1, 

ruo for 1 < 2 < 3, (19) 

nio — 0.7(2; — 3) for z > 3. 



Fitting this model to the quasar luminosity function data in for 2 > 1, Kochanek 
finds that "the best model" has a = 1.07 ± 0.07, b = 0.27 ± 0.07 and = 18.92 ± 0.16 
at B magnitude. The magnitude corrected probability, Pi, for the quasar i with apparent 
magnitude rrii and redshift Zi to get lensed is: 

Pi = T{zi)B{mi, Zi) . (20) 

Selection effects are caused by limitations on dynamic range, limitations on resolution 
and presence of confusing sources such as stars. The survey can only detect lenses with 
magnifications larger than Mj. This sets the lower limit on the magnification. Therefore the 
Ml in the bias function (|T^) gets replaced by Mf{6) (for details, see K93), which is given as 

M^ = Mo(/ + l)/(/-l) , (21) 

with 

/ = ioO-4A™W _ (22) 

The corrected lensing probability and image separation distribution function for a single 
source at redshift zo, are given in K96 

, I .(M)p.(M)B(.^.M,(M).M.) _ 

J B{m,z,Mo,M2) 

and 

Pci - PcA^o) — — — , , , , . — , (24j 

Pi B{m, z, Mo, M2) 

where 

p^{M) = ^-l ^-^dzL. (25) 



t{zs) Jo dzLd{Ae) 

Equation (p^ ) defines the configuration probability. It is the probability that the lensed 
quasar i is lensed with the observed image separation. 

To get selection function corrected probabilities, we divide our sample into two parts — 
the ground based surveys and the HST survey. We use the selection functions as suggested 
in K93. 

In our present calculations we do not consider the extinction effects due to the presence 
of dust in the lensing galaxies. 

Finally, we use the above basic equations to perform the following tests : 
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(i) The sum of the lensing probabihties for the optical QSOs gives the expected number 
of lensed quasars, rii^ = J2 Pi- The summation is over the given quasar sample. We look for 
those values of the parameter for which the adopted optical sample has exactly five lensed 
quasars (that is those values of the parameters for which = 5). 

We start with a two parameter fit. We allow a to vary in the range (0.0 < a < 2.0) and B 
to vary in the range (0.5 < B < 10.0). We observe that for i? > 1, Da becomes independent 
of it. It can be easily checked from equation (|^) that for large values of -B, sinh(a:;) ~ x 
making Da independent of B. From the previous works constraining power law cosmology 



111 , p!2| , the value of -B = 1 is found to be compatible with observations. Incidentally, we 
can estimate the present scale factor of the universe as ao ~ c/i^o, hence we use S = 1 in 
further analysis. 

Fig. |I] shows the predicted number of lensed quasars for the above specified range of a. 
We obtain tzl = 5 for a = 1.06. We further generate 10^ quasar samples (each sample has 
867 quasars) using bootstrap method and find best fit a for each data set to obtain error 
bars on a. We finally obtain a = 1.09 ± 0.3. 

(ii) We also perform maximum likelihood analysis to determine the value of a, for which the 
observed sample becomes the most probable observation. The likelihood function is 

-^=11(1-^:) iiPkPck- (26) 

i=l k=l 

Here Ni^ is the number of multiple-imaged lensed quasars, A^'u is the number of unlensed 
quasars, pj^, the probability of quasar k to get lensed is given by Eq. ^ and p*^, the config- 
uration probability, is given by Eq. |2^. The best fit (-Cmax) occurs for a = 1.13. We see that 
0.85 <a< 1.56 at la (68% confidence level) and 0.65 <a< 2.33 at 2a (95.4% confidence 
level) . 

However, gravitational lensing statistics is susceptible to a number of uncertainties. The 
constraints obtained on the cosmological parameters from the statistics of strong lensing 
may vary after inclusion of the uncertainties in the luminosity function, lensing cross-section 
for galaxies (E/SO) and quasar luminosity function, role of spirals and the dust extinction. 

3.2 Constraints from age estimates of high-z galaxies 

Another observational test which can constraint a is the age measurement of the old high 
redshift galaxies (OHRG)|^ These constraints are more stringent than those obtained 
from globular cluster age measurements p5|, 0, |3^. Here, we consider the galaxy 3C65 at 



z = 1.175 (4 Gyr old) |3, at z = 1.55 (3.5 Gyr old; 53W091) |3|, || and a 4 Gyr old galaxy 
53W069 at z = 1.43 These are the minimum ages of these galaxies as indicated by best 
fitting spectral synthesis models. The age of the universe at a given redshift is greater than 
or at least equal to the age of its oldest objects at that redshift. In power law cosmologies 
the age of the universe increases with increasing a. Hence this test gives lower bound on a. 
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This can be checked if we define the dimensionless ratio: 



t(z) f(a,z) , , 

The tg is the age of an old object and f{a,z) = a/{l + zY^". The error bar on Hq determines 
the extreme value of tg. The lower limit on Hq was recently updated to nearly 10% of accuracy 
by Freedman [Q: Hq = 70±7 km/sec/Mpc; Icr. The constraints from SNe la data also point 
to Hq > 60 km/sec/Mpc For the galaxy 3C65, the lower limit on age (4.0 Gyr) yields: 
0.26 <Hotg< 0.32. Similarly for the galaxy 53W069 we have 0.26 <Hotg< 0.32. For the 
galaxy 53W091 aX z = 1.55 the lower limit on age (3.5 Gyr) gives 0.23 < Hotg < 0.28. 

Fig. H shows the variation of the function f{a,z) with the redshift z for several values of 
a. We see that a should be at least 0.8, in order to allow for these OHRG to exist. 



4 Summary and a combined constraint 

The main results of the present paper along the with constraints obtained from the SNe 
la data [|ll| are summarized in Table 2. The motivation for our work was to establish the 
viability of a linear coasting cosmology a{t) = t. Using gravitational lensing statistics, we find 
that such a coasting is accommodated within la. The age determination of OHRG gives as 
a lower bound a > 0.8 for a power law cosmology a{t) oc t". Moreover, for such a cosmology 
Hoto = a. With updated value of i^o = 70 ±7 km/sec/Mpc and to = 1A±2 Gyr @], it gives 
a = 0.98±0.25. Dev et al (2001) reported that a = 1.0 is consistent with SNe la data (within 
68% confidence level). We find that a = 1.0 is in concordance with the listed observational 
tests. It is interesting to observe that the lensing analysis barely accomodates an Einstein-de 
Sitter universe {a = 2/3) at 2cr. Similarly, the age determination of OHRG also rules out 
a = 2/3. We conclude that the coasting cosmology with strictly linear evolution of scale 
factor, a(t) = t, cannot be ruled out on the basis of these observations. 
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Method 


Reference 


a 


Lensing Statistics 






(i) riL 

f T T 1 1 1 1 T"/^ 1 1 /- X /-"-V *-J 

J ijiKeiinooQ 
Analysis 


This paper 
inis paper 


1.09 ±0.3 

J^-J^ -3-0.3 


OHRG 


This paper 


> 0.8 


SNe la 


Dev et aL(2001) 


1.004 ±0.043 



Table 2: Constraints on a from various cosmological tests. 




0.5 1 1.5 2 2.5 



a 

Figure 1: Predicted number of lensed quasars riL in the adopted optical quasar sample, with 
image separation < 4, Vs power index a. 
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12 3 4 

z 

Figure 2: f{a, z) vs. z for various different values of a; 'a' corresponds 3C65 galaxy [z — 
1.175), 'b' corresponds to the galaxy 53W069 {z — 1.43) and 'c' corresponds to 53W091 
{z = 1.55). 
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